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Abstract 

For mixtures of Bose and Fermi alkali atoms, the Fermion degrees of freedom can be integrated 
out in their finite temperature partition function. The final result can be expanded as a power 
series in the Boson density. Under appropriate conditions, the pairwise interaction between the 
bosons can be changed from positive to negative at a low enough temperature by the Fermion- 
mediated term and the Boson cloud may then collapse. For attractive Fermion-Boson interactions, 
there is also a collapse of the Fermion cloud. This transition is first order (discontinuous) because 
the leading order nonlinear interaction term between the Bosons is third order in Boson density 
with a negative coefficient. We discuss the finite temperature phase diagram of this transition. 
Our result may provide for an explanation of recent experimental observations by Modungo and 
coworkers. 

PACS numbers: PACS: 03.75.Fi, 32.80.Pj 
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There is much recent interest in mixtures of Bose and Fermi alkali atoms in traps 
experimentally. This includes different hyperfine states of the Li6-Li7 system, 

Na23-Li6 and K40-Rb87. The original motivation was to study possible superconductivity 
in the Fermions. In these systems, the Fermions are essentially non-interacting. From the 
Fermi exclusion principle, no two Fermions can be at the same place, the s-wave scatter- 
ing between them cannot occur. As a result, evaporative cooling of the Fermions alone 
is ineffective. To evaporatively cooled the Fermions the Bosons were used as a refrigerant 
experimentally. Much of the current theoretical interest is focused on estimating the pos- 
sible superconducting properties of the Fermion system. Nevertheless, mixtures of Bosons 
and Fermions are interesting in their own right. Modugno and coworkers [f| recently studied 
mixtures of 40 K and 87 Kb with an attractive interaction between the Boson and Fermion and 
found that as the number of Boson is increased there is a "instability value" Nb c at which 
a discontinuous leakage of the Boson and Fermion occurs. They observed a first order 
collapse transition of the Boson and Fermion clouds. This paper explores the origin and the 
nature of this collapse transition. The nature of the transition depends on the higher order 
term of the Boson interaction. Because there is no direct interaction between the Fermions, 
it is possible to formally exactly integrate out their degree of freedom in their finite tem- 
perature partition function. The final result can be expanded as a power series in the 
Boson density. An effective interaction between the Bosons is obtained. We found that the 
pairwise effective interaction can become attractive for a large enough number of Bosons 
and at a sufficiently low temperature and the Boson cloud collapses. This is analogous to 
earlier experiments with fixed attractive Boson interaction [f| where the Boson cloud also 
collapses. In the present case, the effective interaction depends on the atom densities and 
can change during the collapse. In addition, in contrast to the one component case, the 
Boson interaction involve powers of the Boson density higher than the second. The leading 
nonlinear interaction term is proportional to an odd (third) power of the Boson density 
(n 3 B ) with a negative coefficient. As a result the transition is discontinuous and physical 
quantities change discontinuously. Using the experimental parameters, we estimate the 
instability Boson number Nb c for the Boson cloud to collapse. We found Nb c to be ex- 
tremely sensitive to the precise value of the Boson Fermion scattering length cibf- For 
Nf = 2 x 10 4 , Nbc ~ 10 5 when clbf ~ —lQ.7nm, in good agreement with experimental 
results. 
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The Boson-Fermion system may also exhibit a superconducting transition. The Boson 
induced Fermi interaction is attractive and thus the superconducting transition always occurs 
at zero temperature. The superconducting transition temperature for the Fermions is usually 
very low, because it is an exponential function of the dimensionless coupling constant 
a: T sc = T F exp(— 1/a). As a result the superconducting transition is difficult to observe 
experimentally. Thus it is important to know if the same difficulty applies to the collapse 
transition. The renormalized Boson interaction can change from attractive to repulsive at 
a sufficiently high temperature. A phase diagram of N Bc as a function of the temperature 
T is worked out. We found that the collapse transition is much easier to observe than the 
superconducting transition. A possible interpretation of the observation of Modugno and 
coworkers from the above picture is as follows: As a result of the Fermion-Boson attraction, 
the collapse transition will happen to both the Boson and the Fermion. When the density 
is high enough, they disappear from the trap. If the amount of Fermion is smaller than the 
Bosons, they will be gone first. Once the Fermion density is small enough, the Bosons will 
stop collapsing. We now discuss our results in detail. 

We first integrate out the Fermion degrees of freedom in the path integral for the system. 
Our starting point is the functional-integral representation of the grand-canonical partition 
function of the Bose-Fermi mixture. It has the form 

Z = J D[0*]£[0]£[V>*p^ (1) 

0(r, r) describes the Bose component of the mixture, whereas ip{r, r) corresponds to the 
Fermi component. This integral consists of an integration over a complex field 0(r, r), which 
is periodic on the imaginary-time interval [0, ft/3], and over the Grassmann field ip(r,r), 
which is antiperiodic on this interval. The term describing the Bose-gas has the form: 
<p) = ^ dr J dr {</>*(r, r) (h£ - + V B (r) - fi B ) 0(r, r) + f |0(r, r)| 4 } . Because 
the Pauli principle forbids s-wave scattering between fermionic atoms in the same hyper- 
fine state, the Fermi-gas term can be written in the noninteracting form: Sp(ip* = 
Jq 13 dr J dr |V ; *( r ? r ) f^-Jf — + ^M 1- ) — AAf) VK 7 ") r )} • The term describing the interac- 
tion between the two components of the Fermi-Bose mixture is: Si n t(<p* ,<f),ip* — 
9bf Jo d>T J dr\ip(T, r)| 2 |</>(r, r)| 2 , where g B = 4:irh 2 aB/mB and qbf = IntfaBF /mi, mi = 
tubttif / {ttib + mp)- 

The integral over the Fermi fields ip is Gaussian. We can calculate this integral and obtain 
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the partition function of the Fermi system as a functional of the Bose field <f>(r, r). To pro- 
ceed we rewrite the partition function as Z = J D[(J)*}D[(J)}Zf exp j— |- (Sb(4>*, </>))} • Z F is 
the integral involving i/)(t, r): Z F = / D[if)*]D[if)] exp (-± (S F (i()*,il>) + S int (4>\ 0, V))) = 
/D[^]D[^]exp{/ ^rfr/rfr/ ^rfr7rfrV*(r,r)G- 1 (r,r,r , ,r , )^(r , ,r , )} ) where G 1 = 
Gp 1 - SGp^r, r, r', r') = -I (ftJL - |g + V^(r) - 5(r - r')5(r - r');S(r, r, r', r') = 
^|0(r,r)| 2 5(r-r')5(r-r'). 

Using the formula for the gaussian integral over the Grassmann variables 
/ Tin dip n dip n exp j— J2 n ,n' V'nAixVv} = det A = e Tr ^ lnA ^Z F becomes only a function of the 
Bose fields 0: Zp = exp (Tr In (— G _1 ))G contains the additional interaction between the 
Bosons due to the presence of the Fermions. The properties of the Fermions can also be cal- 
culated. For example, the Fermion density can be calculated from the functional derivative 
< ip*(r)i/)(r) >=< bZpjbG~ x {j, r, r, r) >. To illustrate the applicability of this result, we 
consider its contribution to the effective interaction between the Bosons. 

In the absence of the Fermions, the Boson wave function can be represented as a 
wave function <po which is given in the Thomas- Fermi approximation by \<fio\ 2 — {^b — 
Q.bmBOJ B {x 2 + y 2 + \ 2 z 2 ))/ gs- Multiboson interaction can be obtained from an expansion 
of Zr^-G- 1 ) in powers of the deviation of the boson density from |0o| 2; 1 01 2 — <5|0| 2 ~l~ |0o| 2 - 
In this spirit, we write S = E + ST: where S (r, r, r', r') = gBF\4>o(ji r)| 2 5(r — y')5{t — 
t') /ft. The corresponding Green's function is defined as Gq 1 = Gp 1 — S .We have G 1 = 
G 1 -5H = Gq^I - G 5S), TriH-G- 1 )) = Tr(ln(-G *)) - ££=i Tr[(G 5T) n ]/n = 
— Y^=o S n /h.The pair interaction comes from the term S% = Tr(Go5S) 2 /2. The result- 
ing expression can be evaluated in the Thomas-Fermi approximation and assuming that 
the Fermions move faster than the Bose condensates. To illustrate our approach we 
consider the term Si = Tr(Go5S). Let us rewrite Go in the form: Go(r,r,r',r') = 
Z^n-HnWCnWe-^-^/i-ihcu + t n - /i F )ft/3,where u = vr(2 S + l)/ft/3; s = 0,±1,.., 
(-ft 2 V 2 /2m F + V F (r) + ( 7BF |0o(r)| 2 )^(r) = e^r).^ becomes S 1 = TrG S = 
/drdrG o (r,r,r,r)^|0(r,r)| 2 /ft = g BH /(hP) drdr6\4>(r, r)| 2 £ n (r)£(r)/|>7iu; ~ 
(e n - Vf)}= 9BH/(h/3) f^drdr5U(r, r)| 2 X^£„(r)£(r)/ty[e# e *-''*'> + l].In the semiclassical 
Thomas- Fermi approximation one has: 

1 1 f 1 

EU r )C(r) e/3(MF) - 1 =J^KfJ rfp e ^o(P,r)-^) + 1' ( 2 ) 
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where H (p, r) = p 2 /(2m F ) + V F (r) + g BF \<fio(r)\ 2 . We S et 

9BF fW „M2 1 f j_ 1 



5i= ^y j ^i^ r )i (^y rfp e /3(Ho( P , r) -^ +1 ( 3 ) 

Similarly we get 

Bg% F f h P 1 /" e /3(-H"o(p,r)-^t F ) 

S2 = -^mh rfr ^ l0(r ' r)|2)2 (^i dp (i + e^o(P^) )2 - ^ 

where F (p, r) = p 2 /(2m F ) + Vp(r) - #bf|0o| 2 - 

We can write the effective Bosonic Hamiltonian in a power series in x - - |J - 



#gp = / |^-|V0| 2 + K//(r)5|0| 2 + ^(5|0| 2 ) 2 + c^ F (5|0| 2 ) 3 } , (5) 

where 

Ve ff = 0.5D lD g BF J o 1 + e/3(t _ A) , (6) 

„oo e D/2-l g /3(e+V>-A)^ e 
<7e// = <7B - /^D<7 BF ^ (1 ~ e/j(£ _ A)2 • (7) 

A = /if - V F (r) - g BF \Mr)\ 2 , *d = (2m F /h 2 ) D/2 /(2nC D ), C x = 1, C 2 = 2, C 3 = 2vr; 
7D = (2m F /^ 2 ) D / 2 /(27rC^), C[ = 1/2, C 2 = 2, C 3 = 3vr. c = {2m F /% 2 ) D / 2 /{12^C D ). 
For D=3, we find using the low temperature expansion 

Veff = 9BF (l^' 2 + Y^Tj2 K ^ k B T )^j , (8) 

9eff(T) = g B - g 2 BF K 3 fi 1/2 a. (9) 

a = [1 - g(fc B T//2) 2 ]. 

In the same approximation, we get, for D=3, n F (r) = k 3 J* °° 1+ ^1 a) .%endequation and 
at low temperatures we have: 

2 

n F (r) = 2k z 2 /Z + ^L(k B T) 2 . (10) 

We thus expect the Fermion density n F to increase as the Boson density increases. As usual, 
fi F can be determined from the equation N F = J drn F (r). 

Modugno and coworkers^ recently studied mixtures of m K and 87 Rb with an attractive 
interaction between the Boson and Fermion and found a collapse of the Fermion cloud. 
The parameters of the system are the following: a B = 50 A, a BF = -217l|p, the aspect 
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ratio of the trapping potential, A = 0.076. The corresponding Fermi wave vector is k F = 
7.7 x 10 _4 /A, and T F w 360nf^ when the number of Fermions, Nf ~ 2 x 10 4 . 

Eq. (5) looks like a typical Landau free energy expansion of a system with order parameter 
O = 6\<j)\ 2 in a weak external field K//- As we learn from the study of critical phenomena, an 
instability occurs when coefficient of the O 2 term becomes zero. This instability condition 
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of g e ff = imposes a condition on p,. Now g e //(r = 0) = g B — Qbf^^f — QbfI^bI 9b)^^ 
With the Thomas-Fermi approximation [i B = (rriBU B ) 3 ^ 5 (15N B g B /4:ii\) 2 / 5 /2. Substituting 
this expression of \i B into the equation of g e ffij = 0) = 0, we obtain an estimate of of the 
critical Boson number N Bc for the instability of the system: 

N Bc = (87r[(# B 7rm/ / ' g BF m F a BF k F af - l}[i F d B (a B / 'd B ) 3/5 '/ 'gspf (A/15). (11) 



d-B — yfi'/fTtb^b is the transverse harmonic length for the trapping potential of the Bosons, 
a is a finite temperature correction factor defined after eq. (9). The dependence of N Bc on 
a BF is extremely sensitive. In Fig. 1 we show N Bc as a function of a BF for T/T F = 0.3 
in units of nm. N Bc changes by a factor of two when a BF changes by 5 per cent. The 
experimental value of N Bc rs 10 5 suggests that a BF ~ — 16.7nra. 

Even when g e ff < 0, there is still a tunnelling barrier for the collapse to occur. This 
tunnelling barrier diappears when g e ff = — g c where estimates of g c N c has been given pre- 
viously as a function of A. [13] We have solved this equation iteratively, starting with the 
initial guess in eq. (28). We found that the corrections to iV c from this condition is less than 
1 per cent. 

We next address the nature of the phase transition. In previous studies of attractive 
one component Bosons, the free energy is a function of only the first and the second power 
of the Boson density. There are no higher power terms. For the present two component 
system, the effective Boson interaction involve higher power terms of the Boson density. The 
nature of the transition is determined by these nonlinear terms. We found a term of the 
order of cg BF (Sn B ) 3 in the free energy functional. From Landau type arguments, if an odd 
order term with a negative coefficient appears, then at the transition the free energy can be 
further lowered by a finite increase of Sn B ; the transition is first order. This is in agreement 
with the experimental observation of a discontinuous leakage in the number of Bosons and 
Fermions. The usual second order phase transition require a leading nonlinear correction of 
the order of n% with a positive coefficient and no third order term. This is not the case in 
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the present system. 

The experiment was carried out at finite temperatures. Upon increasing the temperature 
the induced interaction between the Bose particles can change from attractive to repulsive. 
From the condition that g e ff(T) = 0, we obtain a phase diagram of Nb c as a function of 
T/Tp. This is shown in Fig. 2. The current system may also exhibit a superconducting 
transition. The Boson induced Fermi interaction is attractive at zero temperature. The 
superconducting transition temperature for the Fermions is usually very low, because it is 
an exponential function of the dimensional coupling constant a: T sc = Tpexp(— 1/a). As 
a result the superconducting transition is difficult to observe experimentally. The collapsing 
transition temperature is much higher and do not have this problem. 

Roth 14] has studied the structure of the mixture numerically at zero temperature with 
a modified Gross-Pitaevskii equation for the bosons which self-consistently includes the 
mean-field interaction generated by the fermionic cloud. Numerical results were presented 
for the instability Boson number Nb c as a function of orf/c^b for different values of as/ds 
for isotropic traps. It is not possible to compare his results to the experiment because 
his results are expressed in terms of a single trap length ds whereas the experimental 
trap is anisotropic and there are two trap lengths that differ by a factor of more than 
three. When dBF/ds changes by a factor of three, Nb c changes by two orders of 
magnitude. 

Viverit and coworkers [l^ have discussed an instability criterion based on a linear stability 
analysis at zero temperature. According to the estimate of ref. (13), the effective coupling for 
the experimental system is repulsive and far away from the instability point g e ff/\gBF\ = 
0.055. Furthermore their instability criterion does not involve Ng. 

If the temperature is higher than the Fermi temperature but below the Bose condensation 
temperature, it is still possible to derive an expression for the effective coupling for the 
Bosons. At high temperatures we have the Boltzman limit (e~ MF /( fcsT ) < 1). In this case 
Eqs. © and ((Zj) take the form: 



Now the Fermi energy is replaced by the temperature. As expected, the renormalization 
effect becomes very small. 




(12) 



(13) 



7 



Eqs. (JHJ) and Q are useful for studying other physical properties of the system. 
For example, for the trapped harmonic potential t4(r) = ^rribUJ 2 (x 2 + y 2 + A 2 ,^ 2 ), the 
interact ion- induced fractional shift to the ideal BEC critical temperature is ATq/Tq — > 
— 1.33(abb/db)(XbNb) 1 ^ 6 where iVj, is the atom number of Bose particles in the trap. Replacing 
the parameters with those obtained from eq. (18) and (19), our result provides for an easy 
way to estimate this change in the presence of Fermions. |l7j ] 

Our result can also be used to calculate the density distribution of the particles of the 
mixture. Using the standard Thomas-Fermi approximation with the renormalized inter- 
action one can obtain the density of the Bose particles. The Fermi density can then be 
obtained from Eq. (fTUj) for T = 0. 

g e ff only depends on the square of the Fermion-Boson interaction gsF- It can become 
negative for a repulsive qbf > 0. In the the Boson density is increased, the Fermion 

density decreases. The system phase separates. After the phase separation is complete g e ff 
becomes positive again. Thus in that case, we can consider a negative g e ff to be a signature 
of the instability towards phase separation. The present calculation provides a unified way 
of looking at these two different phenomena. 

In this paper, we have used a 3D Thomas-Fermi approximation to describe the Fermions. 
This approximation will break down when the anisotropic quantization of the energy levels 
due to the anisotropic harmonic potential becomes important. The larger energy scale 
corresponds to the transverse harmonic frequency and is hu)F,r = 15.16nK. This value is less 
than the experimental temperature (T c = 360nK). We thus expect the anisotropic energy 
quantization effect not to be important in the present case. This quantization effect will be 
important in the investigation of possible superconductivity as the temperature involved is 
much lower. 

In conclusion, we show that because there is no direct interaction between the Fermions, 
the part of the path integral involving the Fermion field can be carried out exactly. A de- 
scription of the Bosons in increasing powers of its density with renormalized couplings is 
obtained. The Fermion-mediated term can rendered the Boson interaction attractive. When 
the Fermion-Boson interaction is attractive, the system collapses. We calculate the insta- 
bility Boson number as a function of temperature and as a function of the Boson-Fermion 
scattering length. The transition temperature is comparable to the Fermi temperature. Be- 
cuase of the structure of the leading nonlinear term in the Boson interaction, the transition 
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is first order. This may provide for an explanation of recent experimental results of Modugno 
and coworkers. 
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FIG. 1: The instability Boson number as a function of the Boson-Fermion scattering length in 
units of nm. 
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FIG. 2: The instability Boson number as a function of temperature normalized by the Fermi 
temperature for a Boson-Fermion scattering length of -16.7 nm for N f = 2 x 10 4 (solid line) and 
N f = 5 x 10 4 (dotted line). 
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